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Abstract. The morphological paths towards equilibrium droplets during the late stages of the dewetting 
process of a liquid film from a liquid substrate is investigated experimentally and theoretically. As liquids, 
short chained polystyrene (PS) and polymethyl-methacrylate (PMMA) are used, which can be considered 
as Newontian liquids well above their glass transition temperatures. Careful imaging of the PS/air interface 
of the droplets during equilibration by in situ scanning force microscopy and the PS/PMMA interface after 
removal of the PS droplets reveal a surprisingly deep penetration of the PS droplets into the PMMA layer. 
Droplets of sufficiently small volumes develop the typical lens shape and were used to extract the ratio 
of the PS/air and PS/PMMA surface tensions and the contact angles by comparison to theoretical exact 
equilibrium solutions of the liquid/liquid system. Using these results in our dynamical thin- film model we 
find that before the droplets reach their equilibrium they undergo several intermediate stages each with a 
well-defined signature in shape. Moreover, the intermediate droplet shapes are independent of the details of 
the initial configuration, while the time scale they are reached depend strongly on the droplet volume. This 
is shown by the numerical solutions of the thin-film model and demonstrated by quantitative comparison 
to experimental results. 

PACS. 68.05.Cf Liquid-liquid interface structure: measurements and simulations - 02.60.Cb Numerical 
simulation; solution of equations 



1 Introduction 

Understanding the dynamics and evolution of patterns in 
micro- and nanoscale two-layer immiscible thin films is of 
fundamental importance for the design of functional sur- 
faces for numerous applications ranging from model for 
tear film of the human eye [Ij to instabilities during the 
fabrication of polymer/polymer solar cells [2j. The past 
decades has seen many investigations of such systems, 
both experimentally and theoretically. But besides of its 
technological relevance and its scientific interest, the in- 
depth understanding of these processes is still superficial 
compared to liquid dewetting from solid substrates and 
quantitative experimental verification is still missing. For 
experimental reasons the morphology and the dynamics of 
liquid/liquid dewetting was typically studied using poly- 
mers as high molecular liquids. For the model two-layer 
system of polystyrene (PS) on polymethyl-methacrylate 
(PMMA), the shape of an underlying liquid (PMMA) and 
the liquid polystyrene (PS) rim profile dewetting from this 
substrate has been studied experimentally in the pioneer- 
ing work of the group of G. Krausch O |4]. As a result, 
they found a characteristic shape and dewetting dynam- 
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ics, depending on the relative viscosity of the two liquids. 
The experimentally observed behavior was claimed to be 
in agreement with another important work by Brochard- 
Wyart et al. [5j. However, this is surprising since the dewet- 
ting velocity strongly depends on film thickness and which 
are not considered in [5j. However, the used polymers in 
[3 , 4J are above the entanglement of the respective chain 
length and viscoelastic properties cannot be ruled out in 
this system. Other groups studied the breakup and the 
hole growth of a liquid-liquid system where the viscosity 
of one of the liquids is much larger than the viscosity of the 
other liquid [6] and in the special case, where the resulting 
dewetting morphologies are coated with a thin layer of the 
underlying liquid [7J, whereas the characteristic shape of 
the liquid-liquid interface was not explored in detail. Sim- 
ilar results on dewetting rim shapes can be found also in 
[6|[8Hl2) for different polymer-polymer dewetting systems. 
However, in none of those experimental studies the pa- 
rameters are worked out well enough to allow for a quan- 
titative comparison with theoretical results. While many 
of the theoretical works have focussed on linear stability 
analysis and numerical simulations of short time and long- 
time evolution [T3^,T!4], Golovin and Fisher [14^, even in 
the presence of surfactant [IGJ, a systematic comparison 
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to experimental results for the simplest possible system 
will be important to build a fundamental understanding 
of this system. 

In this study we present a quantitative comparison of 
experimentally obtained transient and equilibrium drop 
shapes to numerically and analytically drop shapes for a 
model system using short-chained PS/PMMA, that can 
be described via a Stokes model for the multiphase flow 
of two immiscible liquids with capillary boundaries and 
driven by the intermolecular forces between the liquids. 
In our experiments we obtain 3d dynamic and equilib- 
rium droplets shapes by following the complete dewetting 
process. The relevant parameters required for the com- 
parison with theoretical derived dynamic droplet shapes, 
i.e. the surface tension ratio and the contact angle, are ob- 
tained by comparing the experimental equilibrium droplet 
shapes to the exactly know theoretical shapes. These val- 
ues are used as input parameters in our dynamic thin-film 
model to compare the various morphological transitions 
long before reaching equilibrium. Except for a short ini- 
tial phase the morphology was found to be independent 
on the initial start conditions. This independence of the 
transient droplet shapes allows for quantitative compari- 
son of the experimentally determined drop morphologies 
to the theoretically droplet shapes calculated using the 
experimentally determined input parameters. 

2 Methods 

For our liquid/liquid experiments we prepared thin poly- 
styrene (PS) films on top of underlying polymethylmeth- 
acrylate (PMMA) films which are supported by silicon 
wafers: Silicon rectangles of about 2 cm^ were cut from 5"- 
wafers with (100) orientation. The pieces were pre-cleaned 
by a fast C02-stream (snow-jet, Tectra) and sonicated in 
ethanol, acetone and toluene followed by a bath in perox- 
ymonosulfuric acid and a careful rinse with Millipore^^ 
water. On top of the cleaned silicon wafers, PMMA films 
were spun from toluene solution with thicknesses ranging 
from 80 nm to 700 nm. The dewetting PS films cannot be 
spun directly on top of the PMMA with the desired thick- 
ness and were, in a first step, spun from toluene solution 
onto freshly cleaved mica sheets with thicknesses ranging 
from 7nm to 22 nm. In a second step, the glassy PS films 
were transferred onto a water surface (Millipore^^) and 
picked up with the PMMA coated silicon substrates. 

Both polymers were purchased from Polymer Standard 
Service Mainz (PSS-Mainz, Germany) with M^/Mn = 1.03 
and molecular weights of = 9.6 kg/ mol for PS and 
= 9.9 kg/ mol in case PMMA. The molecular weights 
are well below the entanglement lengths of PS and PMMA 
and their melts can be treated as Newtonian liquids [17- 
|2T] at temperatures above their bulk glass transition tem- 
peratures Tg of about Tg,ps = 85 =b 5°C and Tg,pmma = 
115 ± 5°C for PS and PMMA, respectively. The glassy 
sample can be stored in the lab at room temperature and 
the liquid/liquid dewetting process is started when heat- 
ing the sample above the glass transition temperature of 
the polymers. The dewetting experiments were typically 
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Fig. 1. Composition of AFM top and bottom topographies af- 
ter dewetting at T = 140 °C for 45 min. Initial top layer thick- 
ness 80 nm PMMA and bottom layer thickness 20 nm PS layer 
lead to isolated droplets, which are not yet in equilibrium 

conducted at a temperature of T = 140 °C resulting in 
viscosities of /i ps ^ 2/cPas for the dewetting PS and 
Mpmma = 675/cPas for the underlying PMMA film. The 
viscosity values were extracted from [22 , 23j. 

The dewetting process was monitored in situ at 140 °C 
by atomic force microscopy (AFM) in Fastscan Mode^^ 
(Bruker, Germany). To determine additionally the shape 
of the liquid PS/PMMA interface the dewetting process 
is stopped at a desired dewetting stage by quenching the 
sample temperature from T = 140 °C to room tempera- 
ture. At this temperature both polymers are glassy and 
the sample can be easily stored and handled. The glassy 
PS structures are removed by a selective solvent (cyclohex- 
ane, Sigma Aldrich, Germany) and the remaining PMMA 
surface with the shape of the formerly PS/PMMA inter- 
face frozen into it is subsequently imaged by AFM. The 
full three dimensional shape of the dewetting PS struc- 
tures are obtained by combining the surfaces of the related 
topographies as determined by AFM as shown in fig. [T] 



3 Results and discussion 
3.1 Equilibrium droplets 

In this section we aim to characterize the subsequently 
used Si/PMM A/PS/air system by its equilibrium config- 
uration, i.e. liquid PS lenses swimming on liquid PMMA 
supported by silicon. We will first discuss the experimental 
method to measure droplet shapes, give explicit expres- 
sions for equilibrium states and describe how surface ten- 
sion values can be extracted from those and finally com- 
pute contact angles. 

In previous experiments no change of the dewetting 
profile could be detected when cooling and freezing a poly- 
styrene dewetting structure on a solid substrate, e.g. [17\. 
However, in our liquid/liquid system it was found that 
the PS/air interfaces significantly change their shape upon 
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Fig. 2. (left) Cross sections of the same polystyrene droplet 
imaged in situ by AFM while cooling the sample from 
145 °C to 30 °C. The entire duration of this cooling cycle was 
about 70 minutes, (right) Temperature dependent height of 
the droplet as obtained from AFM data with glass transition 
temperatures of PS and PMMA are indicated. The solid line 
shows the expected height change based on thermal expansion 
polymers when assuming a constant width of the droplet. The 
dashed line denotes a similar model calculating assuming just 
the linear expansion coefficients. 



cooling. A PS droplet that raised out of the PMMA by 
e.g. 12 nm at T = 140 °C just raises out of the PMMA by 
1-2 nm at room temperature. And as we aim at a quanti- 
tative comparison with theoretical results it is of utmost 
importance to clarify the origin of the shape change and 
to develop a suitable experimental protocol. 

To observe the shape change of the PS/air interface 
during a temperature change, the PS/air interfaces of sev- 
eral droplets were measured while cooling and heating. As 
a first result, the shape change during heating and cooling 
cycles turned out to be reversible, see fig. [2] for a cooling 
cycle. Considering the viscosities of both polymers as well 
as the experimentally found independence of the droplet 
shape from the cooling/heating rate (tested from millisec- 
onds to hours), a lateral equilibration e.g. according to 
temperature dependent surface tensions is not expected. 
This could be confirmed experimentally as no dependence 
of the in-plane radius of the droplet with temperature 
could be detected within an experimental error of about 
±10 nm. However, the height of the droplet changes signif- 
icantly showing a characteristic step in the shape change 
as function of temperature, see fig. [2] 

Even not knowing the thermal expansion coefficients 
precisely, as slightly different values can be found in litera- 
ture, the height change of the droplet can be explained by 
simple model assumptions. The approximate thermal vol- 
ume expansion coefficients of PS are aps,T<TG ^ ^ • 10~^, 
<^ps,t>Tg ^ 5.5 • 10~^ and for PMMA apMMA,T<TG ^ 
2.7 'lO-"^, apMMA,T>TG = 5.4-10-'^ [24J. The remarkably 
change of expansion coefficient at the glass transition tem- 
perature is characteristic for polymers: At room tempera- 
ture both polymers are glassy whereas the E module of the 
PMMA is larger than that of PS. Accordingly the shape of 
the PS/PMMA interface is expected to remain constant 
during the temperature increase and the entire thermal 
volume expansion of the PS is used to increase the height 
of the droplet. Reaching the glass transition temperature 
of PS, the PS expands almost three times faster with tem- 
perature, i.e. twice as fast as the still glassy PMMA. This 
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Fig. 3. Measured equilibrium drop profile as determined by 
AFM. The height scales in nanometers are shown to the right 
of each panel, (left) bottom profile hi and (right) top profile 

change in thermal expansion coefficient thus explains the 
pronounced step in the droplet height for temperatures 
T > Tg,ps- When reaching the glass transition temper- 
ature of PMMA, the expansion ratios of both polymers 
are again similar. Still assuming that no lateral equilibra- 
tion of the lower droplet shape occurs, the height change 
of droplet with temperature is just slightly reduced as we 
still expect the difference between the linear expansion 
for PMMA and the volume expansion of PS, see solid line 
in fig. |2](right panel). Repeating the same reasoning as- 
suming an isotropic change in droplet volume we have to 
calculate the height change of the droplet based on the 
linear expansion coefficients of both polymers which are 
three times smaller than the corresponding volume ex- 
pansion coefficients. The result of this model calculation 
is plotted as dashed line in fig. [2] (right panel). Comparing 
the result of both model calculations to the experimen- 
tal data we can conclude that the height change of the 
droplets results in fact from the thermal expansion of the 
polymers while the PS/PMMA interface remains about 
constant during this temperature induced shape change. 

Based on this result we validated our experimental pro- 
tocol to obtain quantitative drop shapes: The top side of 
the droplets, i.e. the PS/air interface is imaged at the 
applied dewetting temperature and combined with the 
PMMA/PS interface measured at room temperature af- 
ter removing the PS structure. 

Extraction of surface tension and contact angles 

Using the above validated protocol allows to measure three 
dimensional PS drop profiles on PMMA and to determine 
the contact angles and the ratios of the surface tensions 
from their equilibrium shapes. These values will serve as 
input parameters for the computed drop morphologies and 
could not be found in literature with the needed precision 
[25, 26J. 

The analysis of the equilibrium drop shapes is based on 
the following reasoning: Liquid equilibrium droplets with 
constant surface-tension ctq, deform an underlying liquid 
substrate to minimize the surface energy 

E= [ dr (1) 

a={l,2,3} "^^^ 
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Fig. 4. Sketch of a axisymmetric liquid lens with A, ^2, 
being the PS/PMMA, PS/air, PMM A/air interface respec- 
tively and 7 the triple- junction. Each interface Fi has the 
surface tension cr^. The radius of the lens is a and the con- 
tact angles are Qt^^h (top, bottom). Furthermore for planar 
axisymmetric droplets we have wr^ — —Gr cos Ob — e^sin^b, 



-Br COS Ot + Bz sin Ot, npc. 



for fixed droplet volumes. A classical calculation shows 
that minimizers of E have contact angles determined by 
the well-known Neumann-triangle 



E 

= {1,2,3} 



aaUp^ = 0, 



(2) 



where cTq, > corresponds to the surface tensions of the 
PS/PMMA, PS/air, PMMA/air interface for a = 1,2,3 
and each interface has constant mean curvature H^- 
The normalized vector is tangential to the correspond- 
ing interface F^ and normal to the contact line 7 as indi- 
cated in fig. |4j Introducing the spreading coefficient a = 
(73 — (Ji — (72 one can easily verify that equilibrium droplets 
as in Q only exist if (— cr) > and cri,<j2 > ^(— a). A 
sketch of the equilibrium configuration is shown in fig. [4j 
From a measured equilibrium configuration we can thus 
extract the values for surface tension as follows: If for in- 
stance (T3 is given, then one can determine the other two 
surface tensions from the contact angles by solving the 
linear equation 



cos Of cos ^5 
- sin Of sin 



(3) 



where Ot,Oi, > are measured. For contact angles Ot,Oi, < 
90° a liquid lens has the following equilibrium shape for 
r < a 



h2ix,y) = hoo + 



(4a) 



hi{x,y) = /loo - (^V^r - - V^i" - "'j (4b) 

and hi{x^ y) = h2{x^ y) = hoo for r > a. We use the radial 
coordinate = (x — xq)^ + — ^o)^ and call a the radius 
of the droplet. We used the standard representation of 
interfaces Fi by graphs of functions, i.e. 



Fi = {{x,y,z) e 
^2 = {{x,y,z) e 
^3 = {{x,y,z) e 



z = hi{x,y),0 < hi < /12}, (5a) 
z = h2{x,y),0 < hi < /12}, (5b) 
z = hi{x,y),0 < hi = /12}. (5c) 
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Fig. 5. (left) Top and bottom AFM topography of an equi- 
librium droplet. Two cross sections cut perpendicular through 
the droplet in x-direction (dark symbols) and ^/-direction (light 
symbols) are shown together with the fit using eq. Q (dashed 
line); (right) close-up of the top AFM topography including fit 
function. 



Note that up to an additive constant basically hi{x^y) 
and h2{x^ y) are the morphologies measured by the AFM, 
where h2 is measured at T = 140 °C and hi after the 
removal of the PS by an selective solvent. 

A least-squares fit of Q to the measured AFM to- 
pography data returns the six parameters /loo? 
a, xq and y^. Since both interfaces, i.e. hi and h2 are 
measured independently and the AFM can only measure 
height differences, h^^^ xq, y^ have no absolute value, so 
one might define xq = = and hoo as the values set 
by the preparation of the PMMA layer and as determined 
independently. Also h2 is measured in its liquid state and 
hi after solidification and removal of the PS by a selective 
solvent. Using the values for Ho, and a the contact angles 
can be directly computed via 



= arctan I a/ J 



— arctan 



(6a) 



(6b) 



Inserting this into ([3| gives the surface tensions. 

Fitting spherical caps Q to the top and bottom pro- 
files of several droplets, see ng.|5j we obtain a relationship 
between at, a5 and i^i, i^2, respectively which is shown in 
fig. [gJ For constant contact angles (|6| suggests that the 
relationship between curvature and radius must be linear 
which agrees within the experimental accuracy with the 
experimental data. 

Note that due to the magnitude of systematic and sta- 
tistical errors no evidence for a deviation from a linear 
relation in fig. [6] could be observed. 

From the linear relationship between radius and cur- 
vature shown in fig. |6] and using (|6| we obtain for the 
top angle 6t ~ (1.98 ± 0.07)° and for the bottom angle 
6>b ^ (64 ± 2)° where the error is composed from the 
statistical error in the fit and the systematic error in the 
determination of the droplet radius. Using ([3| we obtain 
the surface tensions of the PS/PMMA interfaces 

0-1 = (0.038 ± 0.002)cr3 = (1.22 ± 0.07) mN/ m, 

and of the PS/air interface 

(72 = (0.984 ± 0.001)cr3 = (31.49 ± 0.03) mN/ m. 
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Fig. 6. (top) Curvature of the top spherical caps and (bottom) 
curvature of bottom spherical caps as a function of droplet ra- 
dius a measured from equilibrium droplets on 400 nm thick 
PMMA substrates (circles) and 700 nm thick PMMA sub- 
strates (crosses) with linear fit (shaded area 95% confidence 
interval) gives H'^ = (1.11 zb 0.02)ab and H'^ = (29 =b l)at 
and corresponding contact angles Ob = (64 it 2)° and Ot = 
(1.98 ±0.07)°. 



based on the PMM A/air surface tension as = 32 mN/ m at 
T = 140 °C taken from [26j. The corresponding spreading 
coefficient is a = (-0.022±0.003)a3 = (-0.7±0.1) mN/ m. 



3.2 Nonequilibrium droplet shapes 

Now we are using the parameters extracted from the equi- 
librium droplet shapes in order to plug them into a thin- 
film model and to investigate the evolution of droplets 
near equilibrium. Therefore we shortly explain how such 
a model is derived from a free boundary problem and 
how the corresponding partial differential equation is dis- 
cretized using the method of finite differences. The ques- 
tion of fidelity of the lubrication approximation for the 
particular combination of surface tensions is beyond the 
scope of this paper and will be investigated in future. Since 
it is quite unrealistic that we can obtain sufficiently pre- 
cise initial data from the experiment and use these for a 
numerical simulation, we first check the dependence of our 
results on a particular choice of initial data. Then we com- 
pare a set of simulations for a fixed initial PMMA layer 
thickness of 80 nm at a specific time with corresponding 
experiments at exactly the same physical time. 

Let us consider the final evolution of PS droplets on 
PMMA as they approach their stationary state. The Stokes 
equation which describes the evolution of the free bound- 
aries A(^), r2{t)^rs{t) by a velocity u or correspondingly 
the evolution of hi(t^ y) and /i2(^, ^, y) has a dissipative 
structure with energy £^ in ([T]) and for Newtonian liquids 
the dissipation 



D 



E 



Mi 



(Vu + Vu^)^(ix dy dz. 



After a standard lubrication approximation, where one as- 
sumes that the interfaces are shallow ||V/ii|p, ||V/i2|P = 
0{e^) for £ <C 1 a formal asymptotic calculation shows 
that hi{t^x^y)^h2{t^x^y) are solutions of a fourth-order, 
degenerate parabolic equation. A standard method to cir- 
cumvent the degeneracy of that equation and to avoid the 



Vy 



dx dy. (7) 



where Vh^ = Vh^{h2 — hi) also has the interpretation of 
an intermolecular potential which encodes van-der-Waals 
interactions and the appropriate Neumann contact angle. 
This leads to the following system of degenerate parabolic 
equations 



{1, 2}, (8a) 



where divergence and gradient are component-wise and 
with mobility matrix 



17,3 
3 "-2 



(8b) 

and /i = /ii//i2 is the ratio of PMMA to PS viscosity. The 
reduced pressures are defined ttj — 5Eh^/Shj. This ex- 
pansion is only valid if £ = y{—oY/o^ 1. All quantities 
above are non-dimensional and scaled using 



[y] 



[t] 



(^2 



for some height scale H = 20 nm which we define in corre- 
spondence to the initial PS layer thickness. Using the val- 
ues of the surface tension and known viscosities mentioned 
in the methods section we obtain e = 0.15, [x] = 133 nm, 
[t] = 2.1 sec, /i = 350 and g\Ig2 = 0.039. Note that in 
[27J some of the authors studied stationary solutions of 
p| and highlighted the connection between £^ in (fTl) and 
E}^^ in ([7|) as an appropriate limiting problem for shallow 
interfaces and as /i* ^ 0. 

For better readability we state the model in non-dimen- 
sional form but all results, e.g. numerical solutions, are 
presented in dimensional form again. The particular po- 
tential we use is 

VhA^) = ^hl/h^ - 3hl/h^ and typically 
we use = 1/64. Note that with such a small it 
is only important that Vh^ assumes its unique minimum 
VhAh^) = -1 and 14^ as > K and < /i* <C 
1. Furthermore our particular Vh^ satisfies Vh^ ^ 00 as 
h ^ 0. With larger instabilities and effects related 
to e.g. spinodal dewetting will become more important. 
Note that for the model reduction to be asymptotically 
correct we need (— cr) = 0{e'^) but cri,<j2,cr3 = 0(1). As 
we showed in the previous section we have < (— a) <C cr2 
but (— cr) ~ (Ji. Nevertheless, we observed, by comparison 
to preliminary calculations based on the full Stokes model 
that the lubrication approximation still yields quite accu- 
rate results. In fact, for many similar thin- film problems, 
see e.g. [28j or [29j lubrication approximation is still useful 
near or even beyond its formal validity. 

In order to solve the axisymmetric system of lubrica- 
tion equations (|8| we use a finite difference scheme with 
spatial and temporal adaptivity. The spatial adaptivity 
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is of second order and heuristically adds points near the 
would-be triple junction, where second derivatives become 
large as /i* ^ 0. The time discretization uses a fully im- 
plicit Euler method and the temporal adaptivity is achieved 
with standard step-size bisection to control the error. With 
the exception of interpolation during a mesh adaptation 
the scheme conserves the mass = 27r hi{t^r)r dr. 
Boundary conditions at r = and r = are Neumann 
conditions for hi and tt^. As initial data for our comparison 
we use 



t= 2 s 



hi{t = 0,r) =m, 
h2{t = 0,r) 



4H + HK 
4H + HK 



h r < To 
r > To 



(9a) 
(9b) 



where the initial disc radii tq and heights h are choosen 
to match the volumes and layer thickness prepared in the 
experiments. 

Now we consider numerical solutions of the thin film 
equation and their evolution towards equilibrium. In par- 
ticular we fist show numerical solutions with two different 
volumes 27r /(/12 - hi)r dr = {1.14, 0.07} pm^. For each 
volume we use two different initial configurations as in (91). 
The strong volume dependence of the approach to equilib- 
rium is demonstrated in fig. [Sj where a synchronization of 
the different solutions can be observed as time progresses 
to infinity. 

For the biggest droplet in the first row of fig. |8] this 
happens only after 45 min, while for the smallest droplets 
in the lower row different initial data are synchronized af- 
ter less than 1 min. In our numerical simulations we have 
also started with asymmetric initial shapes and observe 
a quick approach to symmetrical shapes before the inter- 
mediate stages towards equilibrium have begun, see for 
illustration fig. [7| 



Experiments vs thin-film model 

Using the precise knowledge of the surface tensions of our 
system and the weak dependence on initial data we are 
now able to follow the complete dewetting process to un- 
derstand the stages of equilibration of a droplet by com- 
paring these stages in experiment and simulation. For a 
given measurement at a certain time there is an upper 
bound in droplet size above which droplets are not axisym- 
metric with their shape strongly dependent on initial data. 
In practice this implies that for the used liquid/liquid sys- 
tem of PS(lOfc) and PMMA(10/c) at T = 140 °C and realis- 
tic experimental time scales of a few days we are limited to 
droplets of a few micrometers in radius. For experimental 
reason, as explained in the methods section, it is not pos- 
sible to follow the liquid-liquid interface of a single droplet 
on its way to equilibrium. Thus different droplet volumes 
were imaged at the same dewetting time and compared to 
the related simulations. 

Figure [9] shows this comparison just just after having 
reached the time where the drop shapes are become inde- 
pendent from the initial configuration. The experimental 
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Fig. 7. Numerical solutions of the 2d thin-film equation with 
nonsymmetric initial data and /i* = 1/5 and other parameters 
as before at t = 2 sec, 2 min, 30 min {left) top view [right) 3d 
view showing a quick evolution into a nearly axisymmtric state 
followed by an evolution where /ii ^ 0. 



data compare quite well with the numerical solution of our 
lubrication model for all examined volumes verifying the 
extracted knowledge from the equilibrium droplets. Also 
for the late stages of the experiment shown in fig. 10 this 
proves to be true. 

The remaining differences between simulation and ex- 
periment for larger times are suspected due be caused by 
possible slip effects at the substrate/PMMA interface, or 
intermolecular forces between the substrate and PMMA, 
which we so far have neglected in our model. This will be 
subject of our future investigations. 



4 Conclusion 

In this paper we considered the slow evolution of liquid PS 
droplets on liquid PMMA substrate into equilibrium. We 
extracted surface tensions from equilibrium morphologies, 
which was necessary for two reasons: Firstly, the values 
found in literature depend on temperature and molecu- 
lar weight, whereas the parameter regimes studied in lit- 
erature do not overlap with ours, i.e. T = 140 °C and 
^ 10 kg/ mol for PS and PMMA. Secondly, the sign of 
the spreading coefficient depends sensitively on the precise 
value of the surface tensions and is crucial for quantitative 
comparison. 

The surface tension data obtained by this precise vali- 
dation experiment is then used to compare experimentally 
measured surface/interface profiles with the ones obtained 
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Fig. 8. Numerical simulation of droplet formation and equili- 
bration for two PS volumes and different initial configurations 
at different times, (left) ro = 32L, h = H and ro = 16L, 
h = AH and (right) ro = 8L, h = H and ro = AL, h = AH 
(H = 20nm,L = 133 nm) 



by numerically solving a thin-film equation for liquid two- 
layer films. Since the numerical solutions only depend on 
initial data through the total drop volume, we performed 
quantitative comparisons and found very good agreement 
for droplet shapes at similar times and drop volumes. 
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